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Abstract  —  In  this  paper,  a  modular  Simulink  implementation 
of  an  induction  machine  model  is  described  in  a  step-by-step 
approach.  With  the  modular  system,  each  block  solves  one  of  the 
model  equations;  therefore,  unlike  in  black  box  models,  all  of 
the  machine  parameters  are  accessible  for  control  and 
verification  purposes. 

After  the  implementation,  examples  are  given  with  the  model 
used  in  different  drive  applications,  such  as  open-loop  constant 
V/Hz  control  and  indirect  vector  control.  Finally,  the  use  of  the 
model  as  an  induction  generator  is  demonstrated. 

I.  Introduction 

Usually,  when  an  electrical  machine  is  simulated  in  circuit 
simulators  such  as  PSpice,  its  steady  state  model  is  used;  but 
for  electrical  drive  studies,  the  transient  behavior  is  also 
important.  One  advantage  of  Simulink  over  circuit  simulators 
is  the  ease  in  modeling  the  transients  of  electrical  machines 
and  drives  and  in  including  drive  controls  in  the  simulation. 

As  long  as  the  equations  are  known,  any  drive  or  control 
algorithm  can  be  modeled  in  Simulink.  However,  the 
equations  by  themselves  are  not  always  enough;  some 
experience  with  differential  equation  solving  is  required. 

Simulink  induction  machine  models  are  available  in  the 
literature  [1-3],  but  they  appear  to  be  black  boxes  with  no 
internal  details.  Some  of  them  [1-3]  recommend  using  S- 
functions,  which  are  software  source  codes  for  Simulink 
blocks.  This  technique  does  not  fully  utilize  the  power  and 
ease  of  Simulink  because  S-function  programming 
knowledge  is  required  to  access  the  model  variables.  S- 
functions  run  faster  than  discrete  Simulink  blocks,  but 
Simulink  models  can  be  made  to  run  faster  by  using 
“accelerator”  functions  or  producing  stand-alone  Simulink 
models.  Both  of  these  require  additional  expense  and  can  be 
avoided  if  the  simulation  speed  is  not  critical.  Another 
approach  is  using  the  Simulink  Power  System  Blockset  [4] 
that  can  be  purchased  with  Simulink.  This  blockset  also 
makes  use  of  S-functions  and  is  not  as  easy  to  work  with  as 
the  rest  of  the  Simulink  blocks. 

Reference  [5]  refers  to  an  implementation  approach  similar 
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to  the  one  in  this  paper  but  does  not  give  any  details. 

In  this  paper,  a  modular,  easy-to-understand  Simulink 
induction  motor  model  is  described.  With  the  modular 
system,  each  block  solves  one  of  the  model  equations; 
therefore,  unlike  in  black  box  models,  all  of  the  machine 
parameters  are  accessible  for  control  and  verification 
purposes. 

The  Simulink  induction  machine  model  discussed  in  this 
paper  has  been  featured  in  a  recent  graduate  level  text  book 
[6],  and  it  also  has  been  used  by  different  scholars  in  the 
United  States  [7-9],  Korea,  and  Brazil  [10-12]  in  their 
research. 

n.  Induction  Motor  Model 

The  induction  machine  d-q  or  dynamic  equivalent  circuit  is 
shown  in  Fig.  1.  One  of  the  most  popular  induction  motor 
models  derived  from  this  equivalent  circuit  is  Krause’s  model 
detailed  in  [13].  According  to  his  model,  the  modeling 
equations  in  flux  linkage  form  are  as  follows: 
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(b) 

Fig.  1.  Dynamic  or  d-q  equivalent  circuit  of  an  induction  machine. 


where  d  :  direct  axis, 

q  :  quadrature  axis, 

5  :  stator  variable, 
r :  rotor  variable, 

Fjj  is  the  flux  linkage  (i-q  or  d  and  j=s  or  r), 
vqst  vds :  q  and  c/-axis  stator  voltages, 
vqr)  vdr :  q  and  d-axis  rotor  voltages, 

Fmq,Fmd  '■  <7  and  d  axis  magnetizing  flux  linkages, 
Rr  :  rotor  resistance, 

Rs  :  stator  resistance, 

Xh  :  stator  leakage  reactance  (CO  ,,L/v), 

X/r :  rotor  leakage  reactance  (to  eL/r), 


/  Xh  X,r  J 

iq s,  ids :  q  and  c/-axis  stator  currents, 
iqn  hr  '■  cl  and  d-'dx\s  rotor  currents, 
p  :  number  of  poles, 

J  :  moment  of  inertia, 

Te :  electrical  output  torque, 

7}(or  T/) :  load  torque, 

C0e :  stator  angular  electrical  frequency, 

CO  b  :  motor  angular  electrical  base  frequency, 

CO  r :  rotor  angular  electrical  speed. 

For  a  squirrel  cage  induction  machine,  as  in  the  case  of  this 
paper,  vqr  and  vdl.  in  (3)  and  (4)  are  set  to  zero. 

An  induction  machine  model  can  be  represented  with  five 
differential  equations  as  shown.  To  solve  these  equations, 
they  have  to  be  rearranged  in  the  state-space  form. 


ITT.  Simulink  Implementation 

The  inputs  of  a  squirrel  cage  induction  machine  are  the 
three-phase  voltages,  their  fundamental  frequency,  and  the 
load  torque.  The  outputs,  on  the  other  hand,  are  the  three- 
phase  currents,  the  electrical  torque,  and  the  rotor  speed. 

The  d-q  model  requires  that  all  the  three-phase  variables  be 
transformed  to  the  two -phase  synchronously  rotating  frame. 
Consequently,  the  induction  machine  model  will  have  blocks 
transforming  the  three-phase  voltages  to  the  d-q  frame  and 
the  d-q  currents  back  to  three-phase. 

The  induction  machine  model  implemented  in  this  paper  is 
shown  in  Fig.  2.  It  consists  of  five  major  blocks:  the  o-n 
conversion,  abc-syn  conversion,  syn-abc  conversion,  unit 
vector  calculation,  and  induction  machine  d-q  model  blocks. 
The  following  subsections  will  explain  each  block. 


Fig.  2.  The  complete  induction  machine  Simulink  model. 


A.  o-n  conversion  block 

This  block  is  required  for  an  isolated  neutral  system, 
otherwise,  it  can  be  bypassed.  The  transformation  carried  out 
by  this  block  can  be  represented  as  follows: 
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This  matrix  is  implemented  in  Simulink  by  passing  the  input 
voltages  through  a  Simulink  "Matrix  Gain”  block,  which 
contains  the  transformation  matrix  in  (18). 

B.  Unit  vector  calculation  block 

Unit  vectors  ctwB,,  and  sinQe  are  used  in  vector  rotation 
blocks,  "abc-syn  conversion  block,”  and  "syn-abc  conversion 
block.”  The  angle,  Be,  is  calculated  directly  by  integrating  the 
frequency  of  the  input  three-phase  voltages,  (Df: 

0e=Jco  edt  (19) 

The  unit  vectors  are  obtained  simply  by  taking  the  sine  and 
cosine  of  9,,. 

This  block  is  also  where  the  initial  rotor  position  can  be 
inserted,  if  needed,  by  adding  an  initial  condition  to  the 
Simulink  "Integrator”  block.  Note  that  the  result  of  the 
integration  in  (19)  is  reset  to  zero  each  time  it  reaches  2n 
radians  so  that  the  angle  always  varies  between  0  and  2%. 


block  for  the  current  variables  using  (22)  and  (23),  following 
the  same  implementation  techniques  as  before. 


C=1Vcos0<-+v*sine, 

ids  =  -v„s  sin  6e  +  vds  cos  0f 


(22) 
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E.  Induction  machine  d-q  model  block 

Fig.  3  shows  the  inside  of  the  induction  machine  d-q  block 
where  each  equation  from  the  induction  machine  model  is 
implemented  in  a  different  block. 

First  consider  the  flux  linkage  state  equations  because  flux 
linkages  are  required  to  calculate  all  the  other  variables. 
These  equations  could  be  implemented  using  the  Simulink 
“State-space”  block;  but  to  have  access  to  each  point  of  the 
model,  implementation  using  discrete  blocks  is  preferred. 

Fig.  4  shows  what  is  inside  of  the  block  solving  (1).  All  the 
other  blocks  in  Column  1  are  similar  to  this  block. 

Simulation  Tip  1 :  Do  not  use  derivatives.  Some  signals  will 
have  discontinuities  and/or  ripple  that  would  result  in  spikes 
when  differentiated.  Instead  try  to  use  integrals  and  basic 
arithmetic. 


C.  abc-syn  conversion  block 

To  convert  three-phase  voltages  to  voltages  in  the  two- 
phase  synchronously  rotating  frame,  they  are  first  converted 
to  two -phase  stationary  frame  using  (20)  and  then  from  the 
stationary  frame  to  the  synchronously  rotating  frame  using 
(21). 
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Simulation  Tip  2 :  Beware  of  the  algebraic  loops.  Algebraic 
loops  might  appear  when  there  is  a  feedback  loop  in  a  system, 
i.e.,  when  the  calculation  of  the  present  value  of  a  variable 
requires  its  present  value.  When  Simulink  notices  an 
algebraic  loop,  it  will  try  to  solve  it.  If  it  cannot,  it  will  stop 
the  simulation.  An  algebraic  loop  can  be  broken  by  adding  a 
Simulink  “Memory”  block,  which  delays  its  input  signal  by 
one  sampling  time;  however,  this  might  affect  the  operation 
of  the  system.  If  possible,  avoid  algebraic  loops. 


ks  =';*Jcos0e-v*Jsinee 

[v*  =  Vqs  sin  0  e  +V*COS0C 

where  the  superscript  “s”  refers  to  stationary  frame. 

Equation  (20)  is  implemented  similar  to  ( 1 8)  because  it  is  a 
simple  matrix  transformation.  Equation  (21),  however, 
contains  the  unit  vectors;  therefore,  a  simple  matrix 
transformation  cannot  be  used.  Instead,  vqs  and  vds  are 
calculated  using  basic  Simulink  “Sum”  and  "Product”  blocks. 

D.  syn-abc  conversion  block 

This  block  does  the  opposite  of  the  abc-syn  conversion 


Once  the  flux  linkages  are  calculated,  the  rest  of  the 
equations  can  be  implemented  without  any  difficulty.  The 
blocks  solving  the  rest  of  the  equations  are  also  organized  in 
columns.  The  blocks  in  Column  2  solve  (5)  and  (6). 
Equations  (7-10)  use  the  flux  linkages  to  solve  for  the  stator 
and  rotor  d  and  q  currents.  The  fourth  and  the  last  column 
includes  the  electrical  torque  calculation  from  (11)  and  the 
rotor  speed  calculation  using  the  last  state  equation  ( 12);  the 
implementation  of  which  is  shown  in  Fig.  5.  The  rotor  speed 
information  is  required  for  the  calculation  of  the  rotor  flux 
linkages  in  Column  1 ;  therefore,  it  is  fed  back  to  two  blocks 
in  this  column. 
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Fig.  3.  Induction  machine  dynamic  model  implementation  in  Simulink. 


The  resulting  model  is  modular  and  easy  to  follow.  Any 
variable  can  be  easily  traced  using  the  Simulink  ‘Scope" 
blocks.  The  blocks  in  the  first  two  columns  calculate  the  flux 
linkages,  which  can  be  used  in  vector  control  systems  in  a 
flux  loop.  The  blocks  in  Column  3  calculate  all  the  current 
variables,  which  can  be  used  in  the  current  loops  of  any 
current  control  system  and  to  calculate  the  three-phase 
currents.  The  two  blocks  of  Column  4,  on  the  other  hand, 
calculate  the  torque  and  the  speed  of  the  induction  machine, 
which  again  can  be  used  in  torque  control  or  speed  control 
loops.  These  two  variables  can  also  be  used  to  calculate  the 


Fig.  4.  Implementation  of  (1)  in  Simulink. 


output  power  of  the  machine. 

IV.  Simulation  Results 

A.  Initialization 

To  simulate  the  machine  in  Simulink,  the  Simulink  model 
first  has  to  be  initialized  so  that  it  will  know  all  the  machine 
parameters.  For  this  reason,  an  initialization  file  containing 
all  the  machine  parameters  is  formed.  This  file  assigns  values 
to  the  machine  parameter  variables  in  the  Simulink  model. 
For  example.  Fig.  6  shows  the  initialization  file  for  a  30  kW 
induction  machine.  Before  the  simulation,  this  file  has  to  be 
executed  at  the  Matlab  prompt;  otherwise,  Simulink  will 
display  an  error  message. 

Note  that  this  initialization  file  is  the  only  machine-specific 
part  of  the  Simulink  model  that  needs  to  be  changed  when  the 
machine  being  simulated  is  changed. 

B.  Direct  ac  start-up 

To  test  the  model,  the  induction  machine  with  the 
parameters  given  in  Fig.  6  is  simulated  by  applying  220  V 
three-phase  ac  voltages  at  60  Hz  with  just  an  inertia  load. 


Fig.  5.  Implementation  of  (12)  in  Simulink. 
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Fig.  6.  Induction  machine  model  initialization  file. 


Fig.  7  shows  the  three-phase  currents,  torque,  and  speed  of 
the  machine.  The  machine  accelerates  and  comes  to  steady 
state  at  0. 14  s  with  a  small  slip  because  of  the  inertia  load. 

C.  Open-loop  constant  V/Hz  operation 

Fig.  8  shows  the  implementation  of  open-loop  constant 
V/FIz  control  of  an  induction  machine.  This  figure  has  two 
new  blocks:  command  voltage  generator  and  3-phase  PWM 
inverter  blocks.  The  first  one  generates  the  three-phase 
voltage  commands,  and  it  is  nothing  more  than  a  “syn-abc” 
block,  explained  earlier.  The  second  one  first  compares  the 
reference  voltage,  vref  with  the  command  voltages  to  generate 
PWM  signals  for  each  phase,  then  uses  these  signals  to  drive 
three  Simulink  "Switch”  blocks  switching  between  +V/2  and 
-V/2  (Vd:  dc  link  voltage). 

The  open-loop  constant  V/FIz  operation  is  simulated  for 
1 .2  s  ramping  up  and  down  the  speed  command  and  applying 
step  load  torques.  The  results  are  plotted  in  Fig.  9,  where  the 
response  of  the  drive  to  changes  in  the  speed  command  and 


Fig.  8.  Open-loop  constant  V/Hz  control  Simulink  model. 


load  disturbances  can  be  observed. 

Simulation  Tip  3  \  Do  not  select  a  large  sampling  time.  As  a 
rule  of  thumb,  the  sampling  time  should  be  selected  to  be  no 
larger  than  one -tenth  of  the  smallest  time  constant  in  the 
model.  For  example,  for  a  1  kFIz  switching  frequency,  the 
switching  period  is  1  ms;  then,  a  sampling  time  of  0.1  ms 
would  do  the  job.  Do  not  select  a  smaller  sampling  time  than 
required,  either.  A  smaller  sampling  time  does  not  necessarily 
increase  accuracy  but  the  simulation  time  increases. 

D.  Indirect  vector  control  operation 

By  adding  a  vector  control  block  to  Fig.  8  and  d-q  current 
feedback,  the  induction  machine  can  be  simulated  under 
indirect  vector  control.  The  resulting  block  diagram  is  shown 
in  Fig.  10.  Fig.  11,  on  the  other  hand,  shows  the  vector 
control  block.  Note  that  the  current  feedback  for  the  current 
controllers  comes  directly  from  the  induction  machine  d-q 
model  through  Simulink  “Goto”  and  "From”  blocks.  The 
reason  for  using  these  blocks  is  to  reduce  the  clutter  of  signal 
lines. 

The  results  of  the  vector  control  simulation  are  shown  in 
Fig.  12,  where  the  speed  command  and  load  torque  profiles 
are  similar  to  the  ones  described  earlier.  As  seen  in  this 
figure,  the  speed  tracking  and  the  response  to  the  torque 


Fig.  9.  Simulation  results  —  open-loop  constant  V/Hz  control. 


Fig.  10.  Indirect  vector  control  Simulink  model. 


disturbance  are  excellent. 

E.  Induction  generator  operation 

Recently,  because  of  the  increasing  importance  of 
distributed  energy  resources  such  as  wind  turbines  and 
microturbines,  there  is  renewed  attention  to  induction 
generators.  The  model  described  in  this  paper  can  also  be 
used  to  model  an  induction  generator.  The  only  change, 
compared  with  motoring  operation,  is  that  an  external  mover 
rotates  the  motor  at  a  speed  higher  than  the  synchronous 
speed.  This  change  can  be  implemented  in  Simulink  by  using 
a  negative  load  torque  instead  of  a  positive  one  used  for 
motoring. 

Fig.  13  shows  the  torque -speed  curve  of  the  induction 
machine  in  both  motoring  and  generating  regions.  The 
machine  accelerates  freely  to  almost  the  synchronous  speed 
with  direct  ac  start-up,  and  then  a  large  negative  torque  is 
applied  so  that  the  machine  starts  generating. 

V.  Conclusions 

In  this  paper,  implementation  of  a  modular  Simulink  model 
for  induction  machine  simulation  has  been  introduced.  Unlike 
most  other  induction  machine  model  implementations,  this 
model  gives  the  user  access  to  all  the  internal  variables  for 
getting  an  insight  into  the  machine  operation.  Any  machine 
control  algorithm  can  be  simulated  in  the  Simulink 
environment  with  this  model  without  actually  using 
estimators.  If  need  be,  when  the  estimators  are  developed, 
they  can  be  verified  using  the  signals  in  the  machine  model. 
The  ease  of  implementing  controls  with  this  model  is  also 
demonstrated  with  several  examples. 

Finally,  the  operation  of  the  model  to  simulate  both 
induction  motors  and  generators  has  been  shown  so  that  there 
is  no  need  for  different  models  for  different  applications. 
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Fig.  12.  Simulation  results  —  indirect  vector  control. 
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